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We show that a resummation of infrared logarithms is needed to obtain a sensible 
theoretical description of dijet rates when symmetric cuts are applied to the trans- 
verse energies of both jets. We also present the next-to-leading logarithmic (NLL) 
resummation we carried out for DIS production of two jets selected with the cone 
algorithm. 



, 1 Introduction 



It was observed some time ago '1"2' that next-to-leading order (NLO) QCD calcu- 
lations are not able to describe dijet rates measured experimentally in the whole of 
I the dijet phase space |3l4j . In particular in |2] a study was performed of the cross 

Tj" ■ section o-(A) for producing two jets with transverse energies Eti > Et2, Et2 > i?min 

I and Eti > i?min + A. There it was observed that (t(A) at NLO was finite for A = 

(symmetric cuts), but in the vicinity of that point the slope cr'(A) became unphys- 
i' ically positive and infinite. This pathological behaviour can be attributed to the 

fact that incomplete real- virtual cancellations give rise to large logarithms In Q /A 
Qj I (with Q the hard scale of the process) at all orders in the perturbative expansion for 

• (t'(A). However, an all-order resummation of such logarithms restores the correct 

physical behaviour for the slope, as will be discussed in the following. 

• m 

Ix! 

, 2 Resummation of the jet transverse energy difference 

To understand the physical origin of the divergent behaviour of cr'(A) we first 
observe that this quantity can be related to the differential cross section in the 
transverse energy of the highest Et jet: 

a\A) = --^iEn^E,,i,, + A). (1) 
dEti 

Being the opposite of a physical cross section, the slope has to be always negative, 
as is observed experimentally |3I4| . 

However, if we emit an arbitrary number of soft gluons ki, and define Et.jct = 
I SiGjet ^ recombination scheme that adds transverse momenta vectorially, 

we have, from momentum conservation, 

Eti-Et2 = \ ^ fc„|, (2) 

where is the component of kt parallel to the jet axis (which can be taken as the 
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direction of any of the two hard jets). The definition in eq. ^ imphcitly imposes 
the constraint [ X^i^jets | < A on emitted particle transverse momenta. Such a 
veto on soft radiation is known to give rise to terms a" In™ Q/A, which for A ^ Q 
diverge making <t'(A) become infinite, as was first observed in j2j. 

We then consider the probabihty S{qx) of producing an arbitrary number of 
soft /colhnear (SC) partons such that — X^i^jets ^^li- For A ^ Q, SC emissions 
factorise from the leading order value for the slope <Jq{A) as follows: 

a'(A)~f7^,(A)-S(A), ^{A)= JdqxS{qx)e{A-\qx\). (3) 

Assuming that multiple SC emissions are distributed according to a 'random 
walk' [S] around qx = leads to S(A) ^ A/Q for small A. The resulting slope 
is finite and negative defined as expected. QCD radiation is a particular form of 
random walk, giving rise to the following expression for S(A) (see [S]): 

S(A) = - /"?sin(6A)e-«W, (4) 
I" Jo 

where the positive defined 'radiator' R{b) acts as an effective cutoff on the b integral. 
Its perturbative expansion can be reorganised as follows: 

-R{b)^Lgi{asL)+g2{asL)+asg3{asL) + ... , L^lnb, (5) 

where Lgi{asL) resums leading logarithms (a"i"+^, LL), g2{ctsL) resums next-to- 
leading logarithms {a^L"^, NLL) and so on. The NLL expression for R{b) reads 



R{b) = Rinc{b) + Rsoitib) + Riotib) + Rng{b) . (6) 

^inc (b) collects contributions from radiation that is soft and colhnear to the incom- 
ing parton(s). Its expression is analogous to the well-known Sudakov exponent in 
Drell-Yan transverse momentum distribution \7\. The term i?soft(&) accounts for 
coherence of QCD radiation arising from interference among soft gluons at large an- 
gles. Both i?inc(6) and Rsoft{b) are universal, in the sense that they do not depend 
on the details of the jet algorithm. However, i?soft(^) should be corrected by taking 
into account the fact that not all soft gluons contribute to Efi — Et2, but only those 
outside the jets. The correction is provided by the additional terms i?jet(^) and 
Rng{b), both of which depend on the jet algorithm. i?jct(fo), as well as Rinc{b) and 
Rsoitib), can be computed by simply considering the emission of a single SC gluon 
(including its multiple SC splittings) and exponentiating the result. However, NL 
logarithms arise also when soft gluons inside the jets emit coherently a relatively 
softer gluon outside. Such contributions are typical of non-global observables [H] 
like Efi — Et2, and are embodied in R^glb). Their resummed expression, known 
only in the large- A'^c limit, has been already computed for two outgoing jets in e+e~ 
annihilation, both in the cone algorithm jU] and in the inclusive kt algorithm |10| . 

We now discuss the phenomenological impact of such a resummation for dijet 
production in DIS. The explicit expression for R{b) in this specific case can be found 
in [S]. Here we present numerical results corresponding to Q = 20GeV, xb — 0.01, 
where —Q^ is the virtuality of the photon and xb the usual Bjorken variable. 
Jets are selected with the variant of the cone algorithm proposed in 11 with an 
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opening angle S — 0.3, minimum transverse energy i?min — lOGeV and rapidity 
I77I < 1 in the Breit frame. In fig. Q] we plot the ratio D{A) = a'{A)/aQ{A), with 
the slope cr'(A) computed at NLO in the SC approximation (lower curve) and NLL 
resummed (upper curve). As can be seen, while the NLO curve diverges for A <^Q, 
the resummed curve stays finite, vanishing as A/Q. The numerical analysis of |B] 
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Figure 1. The ratio D(A) = ct'(A)/(t^(A). See ^ for details. 



shows also that the impact of the non-global piece Rng{b) turns out to be negligible 
in the whole range of A values considered. This is consoling in view of the fact that 
its expression is known only in the large- iV^ limit. 

3 Obtaining the total dijet rate 

After a resummation for the slope cr'(A) has been computed, and after matching 
to fixed order calculations, one could obtain the total dijet rate cr(0) by simply 
integrating cr'(A) from the maximum kincmatically allowed value Amax back to 
A = 0. However, resumming the distribution in the transverse energy difference 
is not the only way to compute (7(0). One can envisage the possibility of studying 
other observables, which we will discuss shortly. 

First of all one can observe that in the massless Eq recombination scheme Eti — 
Et2 takes contribution also from the invariant masses of the two outgoing jets. 
This makes the observable global, and therefore within the scope of the automated 
resummation program CAESAR ^]. However, one should be aware of the fact that, 
since a mechanism that keeps En— Et2 small is transverse momentum cancellation 
(see eq. (0)), one should expect a divergence in the resummed distribution provided 
by CAESAR at a critical value Ac (see ^21 for a discussion on this point). 

Another possibility would be to consider any observable V that vanishes in the 
limit of two outgoing partons {En — Et2 is only an example of such an observable). 
One could then resum the differential distribution da(V)/dV and integrate it from 
Vmax to = 0, thus obtaining again the total dijet rate. This has been observed 
already in 2 for the case of the azimuthal angle between the jets, where V = n — cj) 
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was the quantity to be resummed. 



4 Outlook 



What we have done so far |B] shows that we have a clear understanding of the 
physics underlying dijet rates with symmetric Et cuts. As a mandatory test of our 
ideas, we plan to extend our calculation in DIS to the fcf algorithm j so that we 
can compare our predictions with existing data 3 4 . Since the data have very small 
errors, such a comparison would also offer a test of QCD predictions in a three- 
jet environment, complementary to the event-shape measurements proposed for 
instance in |14|| . and, as we believe, less affected by non-perturbative hadronisation 
corrections. Further developments should include the extension of our results to 
other hard processes. In particular a study of dijet rates in photoproduction could 
be exploited to better constrain the gluon distribution g(x) at moderately large x, 
which should be important for the LHC |15| . 
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